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Abstract 

We review the value of the critical exponents f3/i^, and ^/v of ferromagnetic 
Ising model on fractal lattices of Hausdorff dimension between one and three. They 
are obtained by Monte Carlo simulation with the help of Wolff algorithm. The 
results are accurate enough to show that the hyperscaling law df = 2(5 /v + 'y/u 
is satisfied in non-integer dimension. Nevertheless, the discrepancy between the 
simulation results and the e-expansion studies suggests that the strong universality 
should be adapted for the fractal lattices. 

Introduction 

Fractal lattices have self-similar character. Any small part of them is magnified; the 
structure similar to the entire one is observed. They are good candidates to investigate 
the physical properties in non-integer dimensions. At the early of 1980 's, Gefen and 
his coworkers did the pioneer works on phase transitions in fractal networks @, ii §]• 
They studied the ferromagnetic Ising model on fractal lattices by Migdal-Kadanoff bond- 
moving real space renormalization group method and found that the Hausdorff dimension 
should replace the space dimension in the description of the critical phenomena. They 
pointed out that topological effects play an important role at criticality. For example, 
fractals of finite ramification order can be solved exactly and show no phase transition 
at finite temperature; contrarily, fractals of infinite ramification order exhibit a second 
order one. Latter, Bhanot et al. [||, || executed Metropolis Monte Carlo simulations 
and suggested the average number of bonds per spin should be used instead. Bonnier 
et al. H used alternative bond-moving real space renormalization approach and high 
temperature expansion as research tools but could not close the problems. Few years 
ago, two groups |T0| performed Monte Carlo simulations equipped with high efficient 
Wolff algorithm and obtained precisely the critical exponents on Sierpinski carpets. They 
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SP(42,12,1) SP(33,18,1) 

Figure 1: Generators of Sierpinski lattices: SP(42,12) and SP(3^18). SP(33,18,1) embed- 
ded in three dimension space is shown section by section from left to right. 

showed that the hyperscaling law, = 2/3/z/ + 7/z/, is satisfied for two different fractals 
of Hausdorff dimension df equal to 1.8927 and 1.7925. Recently, Hsiao, Monceau, and 
Perreau |]TT| extended their probes to the fractals of dimension between 2 and 3, called 
Sierpinski sponge, and gave the conclusion reinforcing the previous studies. 



Ferromagnetic Ising model on Sierpinski lattices 

We denote the Sierpinski lattice at kth iteration step by the symbol SP(£'^, Nocc, k) 
where £ is the size of the generating cell, d the dimension of embedding space, and N^cc 
the number of occupied sites in the generating cell. The lattice is generated by enlarging 
the {k — l)th iteration one by replacing each occupied site by the whole generating cell. 
Strictly speaking, a Sierpinski lattice is not a true fractal except k tends to infinity. 
When d is equal to 2, Sierpinski lattice is called Sierpinski carpet; and when c? = 3, it 
is named Sierpinski sponge. In fig. 1 we give two generators of fractal lattices SP(4^, 
12) and SP(3^, 18), whose Hausdorff dimension df = logNocc/ ^ogi are equal to 1.792 
and 2.631, respectively. We will put a spin on the center of an occupied site unlike what 
Gefan et al. or Bonnier et al. did to put a spin on each vertex of an occupied site. We 
should notice that putting spins on vertices will cause the ambiguity how to define the 
Hausdorff dimension because the total number of spins is not equal to the total number of 
occupied sites. The Hamiltonian of Ising model on a Sierpinski lattice is defined as usual 
case: H{{si}) = — J2<i,j> ^i^j where < i,j > runs over all the nearest neighbor bonds on 
SP(£'^, Nocc, k) and the coupling constant has been set to 1. We will study fractals with 
infinite ramification order to insure it exhibits a second order phase transition. 



Homogeneity hypothesis and finite size scaling 

The generalized homogeneity of the free energy per spin in integer space dimension d 
was first proposed by Widom JT^ and later by Wilson renormalisation group approach 
||T3| . It states that near the critical point, the free energy per spin behaves as /(t, /i, L) = 
h~^f{W^, hh^^, L/b) under the change of length scale from 1 to b, where t is the reduced 
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temperature {T — Tc)/Tc, h the external magnetic field, and L the size of the system. Due 
to the scale invariant character of fractal lattices, one can extend the above hypothesis 
to fractals by directly replacing the space dimension d by the Hausdorff one df. The 
immediate consequence is that the critical exponents should satisfy the following two 
relations: a + 2/5 + 7 = 2, dfU = 2/3 + 7. According to Fisher's finite size scaling theory 



1^ , for a lattice of finite size L, the specific heat per spin CL{t), the thermal average of 
the absolute value of magenetization per spin niLit), the magnetic susceptibility per spin 
XL(t), and the logarithmic derivative of magenetization 4>L(t) = d\nmLit)/dPB should 
behave like following: 

nLit)=L^-QnitL'/n (1) 

where fl stands for C, m, x, and and equal to a/u, —P/i', 7/1^, and u'^, respectively. 
We can see that the maximum value of physical quantity fl is determined by the function 
Qn{x) with X = tL^^'^ . The critical exponents a/u^ 1 l^i and can be, hence, extracted 
by finding the associated maximum values VL"^"-^ at different lattices sizes. Moreover, if 
the maximum occurs at Tq{L) for a system of size L, the critical temperature for 
a lattice of infinite size can be extrapolated with the help of the equation Tq{L) = 
T,(l + x*^L-^/^). Therefore we can execute Monte Carlo simulations at criticality and 
extract the exponents a/u, P/u, 7/z/, and directly. 



Monte Carlo simulation 

We setup periodic boundary condition on the studied lattices and perform Monte Carlo 
simulation. Wolff algorithm is applied at some temperature for producing spin con- 
figurations. Histogram method is then used to extract all possible information near this 
temperature. We repeat the whole above processes several times to get a sufficient number 
of samples for doing statistical analysis. 



Results and discussions 

In table 1, we list the critical temperatures and the critical exponents z/~^, /?/z/, and 7/z/ 
of the five different fractal lattices. The results of SP(4^,12) and SP(3^,8) are taken from 
the paper of Carmona et al. |]TU1 in which the two lattices are denoted by fractal B and 



fractal A, respectively. The three fractals of dimension between 2 and 3 are reported 
in our recent article The effective dimensions 2/3/z/ + 7/z/ of the five fractals can 

be calculated to be 1.808(40), 1.890(3), 2.640(18), 2.900(23), and 2.955(41), respectively; 
they are consistent with their Hausdorff dimensions. It proves the validity of homogeneity 
hypothesis for fractals. 

The effect of correction-to-scaling can be observed in the studies. It is found to be 
weaker in the case of dimension between 2 and 3 than that between 1 and 2. This can 
be understood by considering the convergence speed of the deviation ratio p{i^, Nocc, k) 
which is defined as the percentage of the difference between the mean number of neighbors 
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SP 


(42, 12)" 


(32,8)" 


(3^18)^ 


(4^56)*' 


(3^26)'' 


df 


1.792 


1.893 


2.631 


2.904 


2.966 




1.078(3) 


1.4813(2) 


2.35090(9) 


3.99893(10) 


4.21701(6) 


1/u 


0.309(8) 


0.59(1) 


1.185(27) 


1.410(36) 


1.503(53) 




0.069(10) 


0.080(1) 


0.3224(14) 


0.493(5) 


0.506(11) 


7/z/ 


1.67(2) 


1.730(1) 


1.995(15) 


1.914 (13) 


1.943(19) 



Table 1: The critical temperatures and exponents of five different fractals. Data repro- 
duced from the article of Carmona et al. [10]. * Results reported in [11]. 



at step k and that at infinite step. p{i'^, Nocc, k) measures the discrepancy between the 
psew do-fractal and the true one and can be shown to be: 



Nr 



f Ns 



where Ns is the number of occupied sites on each surface (boundary) of the generating 
cell and Nj the number of its internal bonds. As k increases, p(£'^, Nocc, k) tends to zero 
faster than the fractals embedded in a two dimensionnal space. 

Physicists believe that the critical properties do not depend on the structure of the 
lattice but only on the dimension of system, the range of couplings, and the symmetry 
of order parameter. This is called the strong universality. According to it, the critical 
exponents of a fractal lattice should fall in the same class like the hypercubic one of 
the same dimension. But comparing with the values calculated by e-expansion [[16 



we 



find that they are not consistent with each other. The simulation results reinforce the 
argument of Gefen et al. that the topology effects of a fractal play an indispensable role in 
determining its critical properties. Strong universality should be revised for self-similar 
lattices. 
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